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Abstract. Wc study the Gross-Pitacvskii equation with a delta function potential, qSo, 
where \q\ is small and analyze the solutions for which the initial condition is a soliton with 
initial velocity vq. We show that up to time {\q\ +Vq)~2 log(l/|(7|) the bulk of the solution 
is a soliton evolving according the classical dynamics of a natural effective Hamiltonian, 
ie + qsech\x))/2. 



1. Introduction 

The Gross-Pitaevskii equation (NLS) with a delta function potential and soliton initial 



idtu + ^d^u — qSo{x)u + u\u\'^ = 

u{x, 0) = e™''^sech(x - ao) , 

offers a surprising wealth of dynamical phenomena. In [11], (and numerically in [12]), the 
authors and J. Marzuola studied the high velocity, Uq 3> 1, case and showed that the 
scattering matrix of the delta potential controls the dynamics. In this paper we describe 
the case of small q. The most interesting dynamics is visible for initial velocities satisfying 
Vq < \q\- The low Vq regime has been studied in the physics literature [3],[H],[2], and the 
behaviour in the intermediate range of g's and v^s, that is between the fully quantum and 
semiclassical cases studied in [11] and in this paper respectively, is still unclear. We state 
the main result here with a slightly more precise version given in Theorem [2] in ^ below. 

Theorem 1. Suppose that in (11. ip we have \q\ ^ 1. Then, on a time interval < t < 
6ivl+\q\yyHogil/\q\), 

(1.2) \Ht, .) - e^-^We^^Wsech(. - a(t))||HMM) < ^^1^1'^'' , 

where a, v, and 7 solve the following system of equations 

(1-3) ^« = ^ ' = -^g5x(sech^)(a) , ^^ = \ + Y~ ~ ^qd^{sech^){a) , 

with initial data (ao,fo,0). 

Compared to numerical results, the theorem gives a remarkably good description of the 
dynamics of a slow soliton interacting with a small delta function potential. For example 
consider Vq = 0, aQ < fixed, and |g| — ^ 0, illustrated in Fig{T] When q < 0, the bulk of 
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q=0.01 , vO=0, a0=-3 
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Figure 1 . The top figure shows the evolution of \u{x, t) | for t>o = 0, Oq = — 3, 
q = —0.01 for < t < 1000. In the bottom figure the dashed curve is the 
computed center of motion, and the continuous curve, the plot of a(t) given by 
(11. 3p . More figures illustrating other cases, some with an even more dramatic 
agreement can be found at ,http : / /math . berkeley . edu/ ~zworski/HZl . pdf ] 



the solution is oscillatory about the origin, with the center moving from ao < to — ao > 0. 
Since 

-V + -qr] (a) = -q?] (ao) , r/(x) = sech(x) , 
the time to complete one cycle of oscillation is 

2 dx 



ao 



|g|^/^A/?7^(x) — ?7^(ao) 
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which is of size comparable to Since the theorem provides an accurate description 

up to time ~ lo'l"^''^ log(l/|Q'|)) it covers many cycles for small enough |g|. When g > the 
soliton is repulsed by the 6 potential and slowly slides to negative infinity with the terminal 
velocity g^/^ - see Fig|31 below. 

The proof of our theorem follows the long tradition of the study of stability of solitons 
which started with the work of M.I. Weinstein [13]. The interaction of solitons with external 
potentials was studied in the stationary semiclassical setting by Floer and A. Weinstein 
[1] and Oh [13], and the first dynamical result belongs to Bronski and Jerrard [1]. The 
semiclassical regime is equivalent to considering slowly varying potentials, 

idtU + ldlu-W{hx)u + u\u\'^ = 0, < /K 1 

0) = e™°^sech(x - ao) , \\W^''^\\^ < C , k<2, 

and that case has been studied in various settings and degrees of generality in [S], [B], [7] 
(see these papers for additional references). The approach of these works was our starting 
point. The results of ^ in the special case of (11.41) give 

(1.5) ||n(t,«) -e^'^We^^Wsech(« -a(t))||j^i(K) < Ch , 0<t< Clog{l/h)/h, 
where 

^a = v + 0{h^), ^v = -hW'{ha) + 0{h^), 
d 1 f 

with initial data (cioi'^o, 0)3 We note that unlike in (11.31) the ordinary differential system 
(11. 6p is not exact - see Fig|2] and the discussion below. 

At first the equations (II. ip and (11.41) appear to be very different: a delta function po- 
tential is very far from being slowly varying. The similarity of (II. 3p and (11.61) is however a 
result of the same underlying structure. As we recall in ^the Gross-Pit aevski equations, 
(II. ip or (II. 4p . are the equation for Hamiltonian flow of 

(1.7) Hy{u)''^^J{\d,u\'-\u\')dx + ^Jv\u\\ V = q6o, V = W{h.), 
with respect to the symplectic form on if^(M, C) (considered as a real Hilbert space): 

(1.8) uj{u,v) = lra J uv , u,v E H\RX) ■ 

When ¥ = 0,7] = sech is a critical value (minimizer) of Hq with prescribed norm: 

(1.9) dSn = 0, S{u) = Ho{u) + ^\\u\\l,. 



I^Strictly speaking the result in [5j describes the dynamics for < t < cq/Zi only. That corresponds to 
smah time dynamics of the potential W. Iterating the full strength of the result of [S] seems to give the 
expected extension to Ehrenfest time \og{l /h)/h [TO] . 
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Figure 2. Comparison of the dynamics of the center of motion of the sohton 
for the Gross-Pitaevskii equation with a slowly varying potential, 

iut = —-Uxx — \u\'^u — sech^{hx)u , h = l/5, /i = l/4, 

and initial condition in (11.11) with vq = 0, cq = —3. The dashed red curve 
shows the solution to Newton's equations used in [1] and [5], the blue curve 
shows the center of the approximate soliton u, and the black dashed curve is 
given by the equations of motion of the effective Hamiltonian 

^ [v^ + sech^{h») * sech^(a)) . 

The improvement of the approximation given by the effective Hamiltonian 
is remarkable even in the case of /i = 1/4 in which we already see radiative 
dissipation in the first cycle. 
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The flow of Hq is tangent to the manifold of solitons, 

M = {e*V^(^-"Vsech(/i(x-a)), a,v,-fER, /i G K+} , 

which of course corresponds to the fact that the solution of (11 .ip with g = and uo{x, 0) = 
e^T+'''o(^-°oVsech(/i(x - ao)), is 

(1.10) u{x, t) = e^7+i«o(x-ao)+^(M2-^.^)^/2^ sech{^{x - ao - vot)) . 
The symplectic form (11.81) restricted to M is 

(1.11) uj\m= fJ'dv A da + vdfi A da + d^ A dfi , 

see §2.4[ The evolution of the parameters {a,v,'y,fi) in the solution u{x,t) follows the 
Hamilton flow of 

rr . _ /^^^ fJ'^ 

J^OlM— — — , 

/ 

with respect to the symplectic form uj\m. 

The systems of equations (II. 3p and (II. 6p are obtained using the following basic idea: if 
a Hamilton flow of H, with initial condition on a symplectic submanifold, M, stays close 
to M, then the flow is close to the Hamilton flow of H\m- 

In our case M is the manifold of solitons and H is given by (11.71) 

(1.12) Hv\m {a,v,j,fi) = ^-^ + -^\V*sech'){fia), 

Z o 2 

and in particular 

HgSo\M= Ho\m sech^ (/ia) , Hw{h,)\M= Ho\m +^fJ.^W{h») * sech^(/i«) . 

The equations (11.31) are simply the equations of the flow of Hgs^^ \m - see §2.51 The equations 
of the flow of Hw{h»)\M are easily seen to imply (II. 6p but some h corrections are built into 
the classical motion. It would be interesting to see if this provides improvement of the 
analysis of [5]. Since our interests lie in the study of various aspects of the delta impurity 
we satisfy ourselves with a numerical experiment which shows that the improvement is 
indeed dramatic - see Figl2J 

In either case, all of this hinges on the proximity of the flow to M and to show that we 
use the Lyapunov function, L{w), introduced in [15] - see §3 Typically, and as is done 
in [5], L{w) is bounded from below so that it controls the norm of w (roughly speaking 
the expression estimated in (II. 2p and (ll.4p ). while {d/dt)L{w) is estimated from above. In 
this paper due to the irregularity of the potential that approach for upper bounds does not 
seem to be applicable but we can estimate L{w) directly, controlling the propagation of 
a, w,7, and /i more precisely. 

The paper is organized as follows. In ^ we recall the Hamiltonian structure of the 
nonlinear flow of (11.11) and describe the manifold of solitons. Its identification with the Lie 
group G = H^kM.^, where H-^ is the Heisenberg group, provides useful notational shortcuts. 
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In [|3]we describe the reparametrized evolution. The starting point there is an application 
of the implicit function theorem and a decomposition of the solution into symplectically 
orthogonal components. That method has a long tradition in soliton stability and we 
learned it from [5]. In ^ we give a self-contained and constructive presentation of well 
known spectral estimates. Weinstein's Lyapunov function is adapted to our problem in ^ 
It is estimated using classical energy. The ODE estimates needed for the iteration of our 
stability argument are given in ^ and a stronger version of Theorem [1] is proved in ^ 

Finally, we make comments on the numerics. The computations of solutions of (11. ip and 
(ll.4p were done using the FORTRAN code described in [T2l §3] and written as part of that 
project by J. Marzuola. Other computations and all the graphics were done using MATLAB. 

Acknowledgments. We would like to thank Jeremy Marzuola for allowing us the use 
of his code for NLS computations, and to Patrick Kessler and Jon Wilkening for generous 
help with various computing issues. The work of the first author was supported in part by 
an NSF postdoctoral fellowship, and that of the second second author by an NSF grant 
DMS-0200732. 



2. The Hamiltonian structure and the manifold of solitons 



In this section we recall the well known facts about the Hamiltonian structure of the 
nonlinear Schrodinger equation. The manifold of solitons is given as an orbit of a semidirect 
product of the Heisenberg group and M+. 



2.1. Symplectic structure. Let be a complex Hilbert space with the inner product 
(•, For W, a totally real subspace oiV {W fMW = {0}), we have V = W + iW ~ W^, 
and we can consider W and V as real Hilbert spaces. 

As a real Hilbert space V is equipped with the natural inner product or metric 

g{X,Y) = Re{X,Y)v, 

and the natural symplectic form 

uj{X,Y) = lm{X,Y)v = giX,tY) . 
In other words g, u, and J, multiplication by 1/i form a compatible triple: 
(2.1) u;{X,Y)=g{JX,Y), g{X,Y) = u;{X,zY) . 

In terms of W"^, we have 



9{X,Y) 



Re A 
ImA 



ReY 
ImY 



UJ 



and 



u;{X,Y)=g{j 



Vl/2 

Re A 
ImA 



Re A 
ImA 



J 



ReY 
ImY 



ReY 
ImY 
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where J is the matrix representing muhiphcation by —i: 

Jl 



J 



-I 



For example, when we consider \/ = C" and W = M", then u is just the standard symplectic 
form. 

In our work, we take V = H^(R, C) C C), and on V we use the inner product. 

The symplectic form u is thus 



(2.2) uj{u,v) = lm J uv , 

and the metric g is 

g{u, v) = Re uv 



Now we consider Hamiltonians and associated Hamiltonian flows. Let if : ^ — > M be a 
function, our Hamiltonian. The associated Hamiltonian vector field is a map : V — > TV, 
which means that for a particular point u E V, we have (S//)^ G TuV. The vector field Eh 
is defined by the relation 

(2.3) ij{v, {Eh)u) = duH{v) , 

where v G T^V, and duH : T^V — > M is defined by 



duH{v) = 

as 



H(u + sv) . 

=0 



In the notation of (12. ip if we use g to define functionals, dHu{v) = g{v, VHu), then {Eh)u 

If we take V = H^(R, C) with the symplectic form fl2.2p . and 

H{u) = J^-\dxu\'-\\u\' 

then we can compute 

duH{v) = Re J {{l/2)dxudxV — |npn?;) 
= Re / {-{l/2)dlu-\u\\)v. 



d2 I 1 2 

,.n — n n 



Thus, in view of (O) and (EJD, 

. 1 



The flow associated to this vector field (Hamiltonian flow) is 
(2.4) u = {^^^)^ = \(\dlu-\u\''u\ . 
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For future reference we state two general lemmas of symplectic geometry: 

Lemma 2.1. Suppose that g : V V is a diffeomorphism such that g*uj = ^{g)uj, where 
/i(^) G C°°(y;M). Then for f eC°^{V,R), 

(2.5) (g-^EfigiP)) = J^f^'f^P^ , P^V. 

Proof. This is a straightforward generalization of Jacobi's theorem which is the case of 
fi{g) = 1. To compute {g~^)^Ef{g{p)), we note 

u,iig~%X, ig~%Ejigip))) = {{g~y u) ,^,^{X , Ef{g{p))) = -l-^uj,^^^(X,EMp))) 

= -j^/fiamix) = -j^^[g*df{pmg-%x) 

and the lemma follows. □ 

Suppose that / G C°°{V;M) and that df{po) = 0. Then the Hessian of / at po, f"{po) ■ 
Tp^y I— > T*^/, is well defined. The Hamiltonian map F : Tp^y — > Tp^ is given by the 
relation 

(2.6) [r(po)X](F) = ^,„(F,FX). 
In this notation we have 

Lemma 2.2. Suppose that N G V is a finite dimensional symplectic submanifold of V , 
and f e C°°(V,M) satisfies 

Ef{p)eTpNcTpV, peN. 
If at Po G N , df{po) = 0, then the Hamiltonian map defined by (12. 6p satisfies 

F{Tp,N)cTp,N. 

Proof. Since is assumed to be finite dimensional we only need to prove the lemma for a 
finite dimensional V (any particular Y G (TpV)-^ can be a value of a vector field in a finite 
dimensional submanifold of V containing N). We can then assume that po = (0,0), and 
that in local coordinates near (0,0), = {(x, ^) | x" = ^" = 0}, x = {x',x") , ^ = {CO, 
•' = (•!, . . . , •jfc), where 2k = dim (see for instance [HI Theorem 21.2.4]). The conditions 
of / mean that 

4"/(x', e', 0, 0) = d^,,f{x', e', 0, 0) = , d/(0, 0) = , 

where we wrote (x, ^) = {x' , C,' , x" , C,") . Hence, the Hessian at (0,0) is given by 

^ ^ ^ ^ ~ L /;,,5"(o,o)_ 
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This means that 

(/" (po)X, Y) = \fXeT,N,Ye (TpN)^ . 

where denotes the symplectic orthogonal. Since the Hamiltonian map, F, is defined by 
{f"{po)X, Y) = u{Y, FY) this proves the lemma. □ 

2.2. Associated symmetries and Noether's theorem. For completeness we comment 
on the Hamiltonian version of Noether's theorem which states that the following three 
statements are equivalent 

EhE = uj{Eh, Ee) = 0, 
E is preserved by the Hamiltonian flow of H, 
H is preserved by the Hamiltonian flow of E. 

For example, consider the mass M = j \u\'^- The associated Hamiltonian vector fleld is 
= iu. We compute 



ct;(SM, Sj^) = — Im J iui{dlu + \u\'^u) = 



The flow associated to H^,/ isu ^ e**M, which is the phase invariance of H and thus solutions 
to dtu = i{d1u + 

Similarly, the time translation, u{x,t) i— > u{x,t + s) gives the conservation of energy, 
H(u), the space translation, u{x,t) u{x + gives the conservation of momentum, 
ImJ UxU. 

2.3. Manifold of solitons as an orbit of a group. For g = (a,v,'j,{x) G M'^ x ]R_|_ we 
deflne the following map 

(2.7) H^3u^ — >g-ueH\ {g ■ u){x) = e'^e'"^''-'''^ fiu{fi{x - a)) . 

This action gives a group structure on M.^ x ]R_|_ and it is easy to check that this transfor- 
mation group is a semidirect product of the Heisenberg group and R+: 

G = i73x]R+, /i ■ (a,w,7) = (-,/it;,7) . 

We recall that the Heisenberg group can be identifled with the group of matrices of the 
form 

1 f 7 

1 a , a, u, 7 e M , 
ij 

and that the semidirect product of H and M+ is deflned by 

(/i, jj) ■ {h', /i') = {h-{fi- h'),fin') , h,h' eH . 
Explicitly, the group law on G is given by 

(a,i;,7,/i)-(a',i;',7',/i') = («",^",7",/i"), 
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where 



v" = v + v'fi, a" = a + — , 7" = 7 + 7' + — , fi" = /i/x' 

Remark. As was pointed to us by Bjorn Poonen, the group acts faithfully on the 4- 
dimensional space spanned by l,f ,0,7 viewed as functions on the group. This can be used 
to see that the group is faithfully represented by the group of matrices of the form 



10 

V fi 

a 1/fi 

7 v/fi 1 



u,a,7 G M, /i G IR+ , 



but we will not use this below. 

The action of G is not symplectic but it is conformally symplectic in the sense that 

(2.8) g*ij = i2{g)uj , g = {h{g) , i^ig)) , ^(^) G M+ , 

as is easily seen from (12. 2p . 

The Lie algebra of G, denoted by g, is generated by ci, 62, 63, 64, 
exp(tei) = (t, 0, 0, 1) , exp(te2) = (0, t, 0, 1) , 
exp(te3) = (0, 0, t, 1) , exp(te4) = (0, 0, 0, e*) , 

and the bracket acts as follows: 

(2.9) [ei, 64] = ei, [eg, 64] = -63, [ei, 62] = -63, [63, •] = , 

so 63 is in the center. The infinitesimal representation obtained from (12.71) is given by 

(2.10) ei = -d^, e2 = ix, 63 = 2, 64 = ■ x . 

It acts, for instance on iS(M) C H^, and by X G 5 we will denote a linear combination of 
the operators ej. 

We have the following standard 

Lemma 2.3. Suppose M. 3 t ^ g{t) is a function and that u G 5(M). Then, in the 
notation of (I2.7p . 

j^g{t)-u = g{t)-{X{t)u), 

where X{t) G g is given by 

(2.11) X{t) = a{t)Kt)ei + 41^2 + m) - a{t)v{t))es + ^^64 , 
where g{t) = {a{t),v{t),-f{t), n{t)). 
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Proof. We differentiate 

g(t) ■ u = exp(i7(t)) exp{—a{t)d^) exp{iv{t)x) exp{9{dx ■ x))u , exp 9{t) = fj,{t) , 
and note that 

dx exp{ivx) = exp{ivx){dx + iv) , 
dx exp{e{d^ ■ x)) = exp{9{dx ■ x))e^dx , 
ixexp(9{dx ■ x)) = exp(9{dx ■ x)){e~^ix) , 
either by direct computation or using (12.91) . The formula (12.111) follows. □ 

The manifold of solitons is an orbit of this group, G -rj, to which Eh, defined in (12.31) . is 
tangent. In view of (12.41) that means that 

for some X E Q. The simplest choice is given by taking X = Az, A G M, so that t] solves a 
nonlinear elliptic equation 

-^V" - + A?7 = . 

This has a solution in if A = yU^/2 > and it then is 7]{x) = /isech(/ix). We will fix 
fi = 1 so that 

ri{x) = sechx . 

Using Lemma [2.11 we can check that G ■ ?7 is the only orbit of G to which is tangent. 
We define the submanifold of solitons, M G Hi, as the orbit of t] under G, 

M = G-r]CHi 

and thus we have the identifications 

(2.12) M = G'-r/~G/Z, T^M = 3 ■ r/ ~ g . 

The quotient corresponds to the Z-action 

(a, V, 7, fi) I— > (a, f , 7 + 2nk, //) , k E Z 

2.4. Symplectic structure on the manifold of solitons. We first compute the sym- 
plectic form uj\m on T^M using the identification (12.121) : 



Since 

rf'{x)dx = 2 , / r]{x)dx'r]{x) = , / dx'r]{x)xr]{x)dx = —1 



we obtain from (I2.10p that 

(2.13) (cjfAf)r,(e2,ei) = 1 , (u;U/)r;(e3, 64) 



□ 
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and all the other (co'fj\/)^(ej, ej)'s vanish. In other words, 

{Lj\M)ri = {dv Ada + d'j A rf/i) (0,0,0,1) = {d{vda + 7c?yu)) (0,0,0,1) • 
To find an expression for uj\nj we use (12. 8p and the following elementary 
Lemma 2.4. If a is a one form on x R+ such that 

0^(0,0,0,1) = (vda + 7rf/i) (0,0,0,1) , g*cr = n{g)a , g , 

then 

a = fivda + 'jdfi . 

We conclude that using the identification (12.1 2p 

(2.14) uj\m= fJ'dv A da + vdfi A da + d'^/ A dfi 

Now let / be a function defined on M, f = /(a,f,7,/i). The associated Hamiltonian 
vectorfield, Sj, is defined by 

uj{-,Ef) = df = fada + f^dv + f^^dfi + f^dj . 

Using (I2.14P we obtain 

(2.15) Ef = ^da + (-^ - ^ V- + + (v- - f,] d, . 
The Hamilton fiow is obtained by solving 

fa Vf^ . fv . r . fvr 

V = , a = — , /i = /-y, 'J = v ff,. 



The restriction of 

H(v^ = 

4 

to M is given by computing by 

/if ^ jj,"^ 

(2.16) f{a,v,-f,fi) = H{g-r]) = — -—, g = {a,v,-f, fi) 

The fiow of (I2.15P for this / describes the evolution of a soliton. 
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2.5. The Gross-Pitaevski Hamiltonian restricted to the manifold of sohtons. We 

now consider the Gross-Pitaevski Hamiltonian for the delta function potential 

(2.17) = ^ I {\dM' - \n\')dx + ^q\um' , 
and its restriction to M = G • r]: 

(2.18) Hq\M= f{a, v,'j,n) = — - + -g/i^sech2(/ia) . 

This is obtained from (12.161) and from calculating 

111 

2?l(^ ■ ^)I(O) = 2^^V(-/^a) = -qfihech^ifia) . 

The flow of {Hg)\M can be read off from (I2.15P : 

fa vf 

V = — = /i^g sech^(/ia) tanh(/ia) 

. /. 

a = — = V 

(2.19) /i 

A = A = 

7 = V— — f\i = -v'^ + -ji^ — g/i sech^(/ia) — -g/i^a sech^(/ia) tanh(/ia) 

This are the same equations as (II. 3p . The evolution of a and v is simply the Hamiltonian 
evolution of (f^ + g/i^sech^(;ua))/2, /i = const. The more mysterious evolution of the phase 
7 is now explained by (12.181) . 

Since /i is constant by the third equation, solving this system reduces to solving the flrst 
two equations. The turning position, Oturiu is given by 

I '3'turr 



sech ^ . 

and Figl3] gives a comparison between atum and the numerically computed turning point of 
the center of the soliton. 

3. Reparametrized evolution 

To see the effective dynamics described in §2.51 we write the solution of (11.11) as 

u{t) = g{t) ■ (ri + wit)) , wit) G H\R, C) , 

where wit) satisfies 

uiwit),Xr]) = 0, VXeg. 

To see that this decomposition is possible, initially for small times, we apply the follow- 
ing consequence of the implicit function theorem and the nondegeneracy of uj\m (see [5l 
Proposition 5.1] for a more general statement): 
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Figure 3. Two plots with q = 0.04 and q = 0.09, respectively, and Oq = 
—10. The blue line is the theoretical prediction of the turning point of the 
soliton, |aturn| = sech~^ [y/y/q), and the red dashed line is the actual soliton 
turning point. For smaller values of q the agreement is outstanding. 



Lemma 3.1. For ^ G/W^ (where the topology onGj'L is given by the identification with 
RxRx X R+) let 

Ut.,5 = {ue Hi : inf ||n - 51 ■ r]\\Hi < 6} . 

geT, 

If ^ ^ ^0 = ^o(S) then for any u G Uy.,5, there exists a unique g{u) G S such that 

uj{g{u)~^ ■u-r],X -7]) = VX G fl . 
Moreover, the map u ^ g{u) is in C^{Uy,,s, S). 
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Proof. We define the transformation 

F : H\R,C)xG — > 0*, [F{u,h)]{X) = uj{h-u-r],X-r]). 

We want to solve F{u, h) = for h = h{u) and by the implicit fuction theorem that follows 
for u near G ■ r} ii for any Qq E G the linear transformation 

is invertible. Clearly we only need to check it for Qq = e, that is that dhF{r], e) : 0—^0*, 
is invertible. But as an element of 0* ®0*, dhF{r], e) = {uj\M)rj, which is nondegenerate. □ 

For § §2.11 and 12.51 we recall that the equation for u (11.11) can be written as 

(3.1) dtU = EH,{u)^ = ^ l{\dM'-\u\')dx + ^q\um'- 
Using Lemma [3. II we define 

(3.2) wit)=g{t)-'u{t)-v, g{t)'^' g{uit)), 

and we want to to derive an equation for w{t). 
By the chain rule and Lemma 12.31 

dMt) = dt{g{t) ■ (r/ + wit))) = git) ■ (r(t)(r7 + w{t)) + d^t)) , 

y (t) a(t)^(t)ei + 4|e2 + m) - a{t)v{t))e, + ^^64 , 
nit) flit) 

git) = (a(t), t;(t), 7(t), yu(t)). Combined with (13. ip this gives 

(3.3) dMt) = -Yit)v - Yit)w + git)^'EH^iigit) ■ (r/ + wit)))) . 
To make this more explicit we apply Lemma [2.11 to see that 

git)-^EH,git) = -j^^a{t)*H, 

(since the action of git) is linear on H^, git)"^ and (5f(t)~^)* are identified). We compute 

g*Hqiu) = ^J (At|c?x(e*^''M(/ix))|^ - ij'^\uiijx)\^)dx + ^g/i^|M(-/ia)p 
1 



(3.4) = - / ifi^\d^^x)\'^ -2vfi'^lmd^uix)uix) + Rev'^fi\uix)\'^ - fi^\uix)\^)dx 



so that 



i/i2g|n(-/ia)p 



1 ^iZ/i^^ _ ..u^ 

-Eg^H^iu) = - ( — —u,,^ + vfiUx - fi^lul^u + —u + nq5i» + na)u ) . 



16 J. HOLMER AND M. ZWORSKI 

For us u{t) = 7] + w{t) and we expand the nonlinear term 

\7] + w\ [T] + w) = T] + 2r] w + T] w + 2\w\ 7] + rfw + \w\ w 

linear quadratic cubic 

Inserting this in (13.31) gives 
Lemma 3.2. If w{t) is given by fl3.2p then 

(3.5) dtw = X{t)w + X{t)r] — ift^Cw + ift^Mw — iqfi5o{x + fia)r] — iqfi5o{x + fia)w . 
where X{t) G g is given by 

def ^ , „.,.N„ '^'„ , / ^ , „.A , „ 



(3.6) X(t) = (-/ia + i;/i) ei 62 + { -i + va - — + — ] 64 

fi \ 2 2 y /i 



and 



£w = —-dlw — 27fw — Tfw + -w , A/'w = 2\w\^rj + 77?/;^ + 



w\ w . 



□ 



We now want to estimate the coefficients of X{t) in (13.50 using the symplectic orthogo- 
nahty of Yrj, Y & q and w. For that we define 

P : 5'(M,C) — > Q 

as the unique hnear map satisfying 

uj{u- P{u)r],YT]) =0 VF G0. 

We will need the following 

Lemma 3.3. Let || • || be a norm on g obtained by using the standard norm in the basis 
given by (I2.10p . Then for w G H^, and F G g, 

\\P{Yw)\\<C\\Y\\\\w\\l2, 

\\p{^^^u)\\ < c\\w\\l2 (1 + Iklllilklli) , 



|P((25o(a^ — 2:0) w) II < C||w|||^i ||w" ^ 



1 J. 

/fi|l'"^llL2 ' 



with the constant independent of xq. 



SLOW SOLITON INTERACTION WITH DELTA IMPURITIES 17 

Proof. We start with an explicit expression for P which follows from (12.131) : 

4 



Pi{u) = —Lu{u,e2'ri) = Re y u{x)xri{x)dx , 
^'^''^'^ P2{u) = uj{u, eiT]) = — Im y u{x)dx'r]{x)dx , 

Psi^) =a;(u,e4?7) = Im J u{x)dx{xr]{x))dx . 
Pi{u) = —uj{u,e'irf) = Re / u{x)ri{x)dx . 



We now recall that ||^^||^cx)(i[j) < C'||^^||L2(R)||^^||/i'i(R) and the estimates follow. □ 
Since Pw = Pwt = 0, (13.51) gives 

(3 X,{t) X{t) - g/iP(z5o(« + afi)v) 

= -P{X{t)w) + ^i^PiiCw) - fi^PiiAfw) + qiJ.P{i6o{x + ^ia)w) . 

Since C is the Hessian of given in (11.91) . at the critical point 77, and 'Eg is tangent to M, 
Lemma 12.21 (or a direct computation) shows that 

P{iCw) = , 

and hence that term can be dropped from the right hand side. We can then use Lemma 
13.31 to obtain 

Proposition 3.4. Suppose that w{t) is given in Lemma \3.S\ and that Xi{t) is given by 
1KB. Then 



□ 



(3.9) ||Xi(t)|| < Cq\\w\\H^ + C{\\w\\l, + \\w\\j,,) . 



Finally we interpret the coefficients of Xi{t). First we use (13.71) to see that 



P(i(5o(» + a/i)?7) = -d^{ri )(a/i)e2 + yq (a/i) + -afid^{r] ){aii) ] 63 
Then we combine this with (13. 6p and (13.81) to obtain 

-^iW = + vii) ei + [ -]-qiidx{rf){aii) - - ) 62 



, 2^' ^ ' ' ' u 

^ 1 V- , 

—liqrfiaiJi) qan'^dx('r]'^)(a(x) — 7 + fd 1 les 64 . 

/i 
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We now see that 
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Xi{t) = equations IHIM hold. 



4. Spectral estimates 

In this section we will recall the now standard estimates on the operator C which arises 
as Hessian of £^ at 77: 

Cw = — dlw — 2r]'^w — rfw H — w , 

or 

^ TL I 1 [Re w] r 1 , 1 \ 2 1 

Cw = ^ J T , ^± = — di - (2 ± iW + - . 

In our special case we can be more precise than in the general case (see [I5], and also [5l 
Appendix D]). The self-adjoint operators L± belong the class of Schrodinger operators with 
Poschl- Teller potentials and their spectra can be explicitly computed using hypergeometric 
functions - see for instance [9|, Appendix]. This gives 

cr(L_) = {0} U [1/2, 00) , a{L+) = {0, -3/2} U [1/2, 00) . 

The eigenfuctions can computed by the same method but a straightforward verification is 
sufficient to see that 

L_r/ = 0, L+(9.r/) = 0, L+{r]') = —r^' . 

We now have 
Proposition 4.1. Suppose that for every X G g 

uj{w,X -ri) = 0, weH\R,C)- 
Then, with {w,v) =^ Re / wv on H^{W,C) (considered as a real Hilbert space), 

(4.1) (£^,^/;)>po||^||i2, Po = ^^Y^^p^^ 0.2058. 

We need the following elementary 

Lemma 4.2. Let V be a real vector space with an inner product (•,•), and let L be a 

symmetric operator on V. Suppose that for vq,Vi G V, \\vj\\ = 1, we have 

Lvq = -CqVo , Co > , {Vo, Vi)^ = 02 , 

(4.2) 

{w,vo) = =^ {Lw , w) > ci\\w\\ , ci > . 

Then 

(4.3) {v,Vi)=0 {Lv,v) > C3\\v\\^ , C3 = CiC2-Co{l-C2). 
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Proof. For reader's convenience we present the straightforward argument in which we can 
assume that < C2 < 1. For v ^ V we write v = avo + w, {vq.w) = 0. The condition 
{v,vi) = gives 

(4.4) = -{w,vif = -{w,vi - clvo)^ < 

Hence 

{Lv, v) > CiWwW^ — Cott^ 

>c,6\\wr+(^c,{l-6)^^^-co^ a' 

= (C1C2 - Co(l - C2))||w|P , 

if we choose Ci6 = (ciC2 — co(l — C2)). □ 



Proof of Proposition \4-l\ The assumption means that 

{ir] \ 
.^^ } dx = . 
1X7] 
dxixr]) J 

Working with real and imaginary parts the proof reduces to lower bounds on L±: 

{v,r]) = {v,XT]) = =^ {Lj^v,v) > po\\v\\l2 , 
,77) = {v, {xd^ + 
where now v G iJ^(M; M). Noting that 



(4.5) „ 

{v, d^T]) = {v, {xd^ + l)r/) = ^ {L_v, v) > Pq\\v\\^2 , 



we can apply Lemma 14.21 in the following three cases: 



2 



2 ' 2 ^ 32 

V = {r]Y n H\R,R) , vo = ^dxrj, = ^^xtj , L = L, 

V2 vr 

2 ri 2 ^ 2 ^ 

Cq = U , ''1 ~ 2 ' '"^ ~ ' 

1 2 /2 



V2" ' V^+TT 



Cq — , — , C, 



2' " 12 + 7r2 



20 J. HOLMER AND M. ZWORSKI 

Here we used 



2 

6" 



f sech^(x)(ix = 2 , f sech'^{x)dx = - , f sech'^(x)dx = — , f x^sech^(x)(ix = — , 
Jr Jr 3 Jig 2 Jjg 

I tanh^(x)sech^(a;)(ix = - , f {dx{xsech.{x))Ydx = :^(12 + tt^) . 
Jr 3 Jr lo 

It follows that we can take 

'Stt^ 3 9 9 \ 9 



... /3;]-2 
Po = ^.min^(c{4 - c^o(l " 4)) = "^^^ ( 



16 2' 27r2' 2(12 + 7r2); 2(12 + vr^) ' 
completing the proof. □ 
Proposition 14.11 gives a slightly stronger statement: 
{Cw, w) > {1 — 6){Cw, w) + 6po\\w\\'j^2 



>(l-5) (^l\\dM\'-l\\wr^+Spo\\w\ 
>T^^\\dxwf -0.0760\\d,w\W S 



L2 



5 + 2po n . n , 5 + 2po 

In addition, 

(4.6) {Cw,w) > ~ 0.0555||u;||^i . 

7 + 2po 

Remark. The smallness of these constants gives a possible explanation of the size of g's 
for which the asymptotic result agrees with numerical simulations. The implicit constants 
in ^ are closely related to the constants above. 

5. Estimates on the Lyapunov function 

Suppose u = u{x, t) solves dlll]) witlH I g| ^ 1 and initial data 

(5.1) «o(x) = e"^«r/(x-ao), |t'o| < 1 

Let T > be the maximal time such that on [0,T], the smallness condition 5 < 5q m 
Lemma 13.11 is met. From Lemma 13.11 obtain the parameters /i = /i(t), 7 = 7(t), 
V = v{t), a = a{t) satisfying the symplectic orthogonality conditions stated there. Let 
u = u{x, t) be defined by 

(5.2) u{x, t) = git) ■ u{x, t) = e*^e^^>u(p(x - a), t) , 
and let 

w{x, t) = u{x, t) — r]{x) . 



^The symbol <C 1 means smaller than an absolute positive constant, i.e. one independent of all parameters 
in this problem. 
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The Lyapunov function of [15] and [5] is given by 

(5.3) L{w)=8{t] + w)-8{t]). 

The lower bound on L{w) follows from the spectral estimates of §11 and in particular 
from (14.61) . For the upper bound we will use the conservation of Hq{u) and its relation to 
8{r] + w). 

For future reference we state the following crucial consequence of the orthogonality con- 
ditions on 1/7, and in particular of the condition that Im j irjw = Re / wrj = 0: 

Lemma 5.1. Suppose that for every X G g 

uj{w,X -ri) = 0, weH\R,C). 

Then 

(5.4) - 2 
Proof. We first compute 



\ri + w\ 
Lservati 

lemma Re(w,?7) = and hence 



where we used the conservation of the norm. As noted before the statement of the 



\v + w\\l,=2 + \\w\\l,, 



from which the conclusion follows. 



□ 



As a consequence, we can dispense with /i in the estimates of Proposition 13.41 and we 
reformulate it as 

Proposition 5.2. Suppose 1 — /i <^ 1 and |g| < 1. Then 

\v -a\ + \v + qd^rf{a)/2\ + | - qrj^{a) - qad.^rf{a)/2 - 7 + + 1/2| 



Proof. We use (15. 4p in (13.91) . For example, 



-qd,^rj^{a) + v 



< 2 



-^d^r^\a) + - 
2 n jj, 

l;qfid^T]'^{afi) + - 
2 /i 

1 n 2/ \ 

-qfid.^r] (a/i) + - 



c|g| |1 — /i| 



We also use the estimate for \v — a\ to replace va by in the equation for 7. 



□ 
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We adopt the following notational convention: denote the initial (time t = 0) configu- 
ration of the system by 0-subscripts - uq = m(0), wq = w{0), and ao = a(0), vq = f (0), 
/Xo = /^(O), 7o = 7(0). Similarly, denote the configuration of the system at some fixed time 
ti by i-subscripts. Finally, the configuration at any arbitrary time t we denote without 
subscripts - w = w(t), u = u(t) and a = a(t), v = v(t), fj, = fj,{t), 7 = 7(t). 

With this notation we now state 

Lemma 5.3. Suppose fiQ = 1 and Wo = (equivalently, suppose (15. ip holds), and suppose 
that T > is the maximal time for which the smallness condition in Lemma \3.1\ holds. 
Suppose that for an interval of time C [0,T] , the following conditions hold 

< 1 — /.i ^ 1, max \v{s)\ <^ 1, < 1 , 

(5.5) — 

- til < 1, \ti+i-ti\ max |t;(s)|<l. 

Then there is an absolute constant c* > 1 such that 

sup ||u^(s)|||'i < c^lltfilll-i + c*|g|^ . 

We remark that the inequality, < 1 — yU, in (15. 5p is not an assumption but follows from 
Lemma 15.11 

The main result of this section is the following consequence of this: 

Proposition 5.4. Suppose /io = 1 and Wq = 0, and suppose that T > is the maximal 
time for which the smallness condition in Lemma 13. 1\ holds. Let 

^^^ S\og{l/\q\) ^ 

log 

and suppose there is a partition of the time axis 

= to < ti < ■ ■ ■ < t„ < T 
such that on each subinterval [tj,tj+i], (15.51) in Lemma \5. 3\ holds. Then, 

sup \\w{s)\\l,<\q\''~' 

0<s<t„ 

Proof of Lemma \5.3\ We start by noting that in the argument that follows, we will not use 
any information about w or the parameters /i, 7, a, and v for times < t < t^; only that 
/io = 1 and wq = 0. 

We will conveniently reexpress L{w) given by (15. 3p using the conserved Hamiltonian and 
mass. Since u = g ■ u, we recall (13. 4p to obtain: 

(5.6) Hg{u) = g* Hg{u) = ^1^""^ j \u\^ + ^ft^vlm j d^uu +'^^^ j \d^u\^ 

If 1 
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The expression for the mass, M{u) = J \u\'^, becomes M{u) = fJ' J \u\'^- Using this and 
(15. 6p . we obtain 

Now substitute u = r] + w and use the orthogonahty condition Im f wdxTj = to obtain 

1 1 

(5.7) £iv + w)= j,H,{u) + —M{u) - (^M(n) + 



Im / wdxW ?7(— /ia) Re ty(— /ia) |w;(— /ia)|' 

2n J jj 2/x 



Note that the classical energy term (with the /i terms dropped) 

E{u) \v'M{u) + ^qrjiaf , 
has appeared in this expression. Evaluate (I5.7p at t = ti to obtain 

(5.8) £{v + w,) = \h^{u) + -^M(u) - (f^M(«) + -^//(-/..a.) 

- -^Im Wid^Wi - —r]{-fiiai) Rew{-fiiai) - ^\w{-fiiai) 

By taking the difference of the right hand sides of (15.71) and (15.81) . we obtain 
£{7] + w)- £{?]) 



^?-^f)^'"'' + 4t-^ 

2 ^ . , „,2 



Im f wd^w — —ri(—jj,a)Rew(—jja) — —\w(—j2a)\'^ 

2fi J fi 2/i 

+ -^Im / Wid^Wi^ r]{-Hiai)ReWi{-iJiai) + — \wi{-Hiai)\'^ 

J fii Zfii 

+ {£{r] + w,)-£{v)) 
= I + II + III + IV + V 

where each line has been labeled by a Roman numeral. From the spectral estimate Propo- 
sition |1]T] (see (I4.6P ). we have 

(5.10) cilltfll^i - \\w\\]^i - 11^1 < ^{V + 'w) - £{vi) 
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We next estimate the right-hand side of (15.91) . hne by hne. For ti <t < tj+i, let 

e{tf' = sup 

ti<S<t 

Estimate of the 1st line of (15. 9p . By the assumption wq = and /iq = 1, we have 

(5.11) M{u) = M{r])=2 

and 



(5.12) 



By substituting (I5.1ip and (15.121) into Term I, we obtain I = !„ + If,, where 



and 



I. 



1 n 



1 f 1 



6V/x3 2V/i fii 



\q\ 



Inserting (15.41) in Term I^, gives 
1/1 1 



6 V/i /i, 

6 V/i /ij 
1 „ 



1 - 



1^ 1„ 1,2 1„ „2 



L2 



1 



W\\l2 



(1 



1 



and thus 

For Term I;,, we have 



1 

6' 



Hal < ^(Mh + \\Wi\\l,f 



1 1 



1 \ /I 



fi Hi/ \fi'^ Hfii fif 



\q\ 



and thus 



Collecting these estimates, we obtain 

(5.13) |I|<6^ + 2(^;o' + |g|)6^ 

Remark: This direct calculation is in fact the consequence of dSr/ = 0. We are using 
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which follows from 

Estimate of the 2nd line of f l5.9p (classical energies). We compute 



and thus by Proposition 15.21 

+ l^'(a)) I < c{\\w\\l. + Iglll^ll^. + \\w\\U{\v\ + |g|) 
By the fundamental theorem of calculus, 

As in the proof of Proposition 15. 21 we can install /i's in this expression using (15.41) to obtain 

(5.14) |II| <c{e^ + \q\e){\t-t,\ max \v{s)\ + \q\\t - ti\ + v"^ + \q\) 

ti<S<t 

Estimate of the 3rd and 4th lines of (15.91) . By the Cauchy-Schwarz inequality and the 
Sobolev embedding theorem, 

|III| < \v\\\w\\l^ + \q\\\w\\m + \q\\Mli 

Similarly, 

|IV| < \v\\\w,fH^ + \q\\\wi\\Hi + \q\\\wi\\li 

and thus 

(5.15) |ni| + |IV| < 2{\v\ + \q\)e^ + 2\q\e 
Estimate of the 5th line of (15. 9p . By definition of we have 



(5.16) S{r] + Wi) 



^ j IdxTj + d.:^Wi\^ j \ri + Wi\^ + ^ j \ri + Wi 



Substitute into (I5.16P the three expansions: 

\d^ri + d^Wi\^ = Id^r]]^ + 2Red^r]d^Wi + \dcWi\^ 

\r] + Wi\^ = r]^ + 4:Rer]''^Wi + 2?7^(2(Re ty^)^ + 1^;^^) + Ar]{ReWi)\wi\'^ + \wi\^ 
\ri + ujjp = r/^ + 2?7 Re Wi + \ wi\'^ 

and observe that the linear terms cancel since rj solves —\ri + ^ri" + rf = 0. Thus, we obtain 
the estimate 

(5.17) |V| < 8\\wi\\l-, + AWwiWir + ||u;,||^i < 10||W7,||^1 
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This completes the hne-by-hne estimation of the right-hand side of (15.91) . By combining 
flSlTOj) . and the estimates (l5:T3|) . (l5:TD . fl535D . (l5Tn) for the right-hand side of ([ESD, we 
obtain 

Cie^ < + \e'^ + cfe^ + |g|e + e^)(|t - tA max \v(s)\ + \q\\t - tA + v'^ + \q\) 
4 ti<s<t 

+ [e' + 2{vl + \q\)e^] + [2{\v\ + \q\)e^ + 2\q\e] + 10\\w,\\l^ 

By hypothesis, every term on the right side has a small coefficient, and thus can be 
absorbed on the left side. Therefore, we obtain 



By applying the Peter-Paul inequality \q\e < |c|gp + we obtain the desired estimate. □ 



Proof of Proposition \5.4\ Now let 

e'^{t) = sup 

0<s<t 

On the first interval [0,ti], we apply Lemma ESI with i = 0, and since Wq = 0, we obtain 

e(ti)' < c,\q\^ 

On the second interval [^1,^2], we apply Lemma E3] with i = 1, and since < c^,|gp, 

we obtain 

e{t2f < ic, + cl)\q\^ 
We continue, and after the n applications, we obtain 

eiQ' < c. (^gcij \q\' = c. (1^) \q\' < c^^lgp 

Since we want c^'^^q'^ < Iql"^'^, we require 

„^ «og(l/k|) 



log 

□ 



6. ODE ANALYSIS 

The assumptions of Lemma 15.31 involve estimates on v{s). To control these we use 
Proposition 15.21 and ODE estimates which we present in this section. 

Lemma 6.1. Suppose q is a constant, \q\ <^ 1, and a = a{t), v = v{t), ei = ei(t), 62 = ^2{t) 
are real-valued functions. Suppose f : W ^ is a C"^ mapping such that \ f\ and \ f'\ 
are uniformly bounded. Suppose that on [0,T], 

( d = V + ei a(0) = Oq 

^^■^^ \^ = g/(a) + e2' viO) = vo 
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Let a = d{t) and v = v{t) be the real-valued functions satisfying the exact equations 

a = V a(0) = ao 

^ = qfia) ' ^(0) = Vo 



en 



(6.2) 



with the same initial data. Suppose that on [0,T], we have \ej\ < \q\ for j = 1,2. Th 
provided T < ^Igl"-*^/^ log(l/|g|), we have on [0,T] the estimates 

|a-a| < |g|^"^'^log(l/|g|), \v - v\ < Iglt"^"^ log(l/|g|) 

Before proceeding to the proof, we recall some basic tools. 

Gronwall estimate. Suppose b = b(t) and w = w{t) are real- valued functions, g is a 
constant, and (b,w) satisfy the differential inequality: 

\b\ < \w\ b{0) = bo 

\w\ < \q\\b\ ' w{0) = Wo 

Let x{t) = \q\^/^b{\q\-^^H), y{t) = w{\q\~^^H). Then 

\x\<\y\ x{0) = xo = \q\^^X 

\y\ < \x\ ' y(0) =yo = Wo 

Let z{t) = x"^ + y"^. Then \z\ = \2xx + 2yy\ < 2|a;| |y| + 2|a:;| |?/| < 2{x^ + y^) = 2z, and hence 
z(t) < 2(0)e2*. Thus 

\x{t)\ < V2max(|xo|, |?/o|)exp(t) 

\y{t)\ < V2max{\xo\,\yo\)exp{t) 
Converting from {x,y) back to {b,w), we obtain the Gronwall estimate 



exp(|^ 


r|V2t) 






L/2 



(6.3) 

\w{t)\ < V2max(|g|i/2|^o|,ko|)exp(|g|i/2t) 

Duhamel's formula. For a two- vector function X(t) : M — R^, a two- vector Xq G M^, 
and a 2 X 2 matrix function A{t) : M — * (2 x 2 matrices), let X{t) = S{t,t')Xo denote 
the solution to the ODE system X{t) = A{t)X{t) with X{t') = Xq. In other words, 
-^S{t,t')Xo = A(t)S{t,t')Xo and S{t',t')Xo = Xq. Then, for a given two-vector function 
-F(t) : M ^ M^, the solution to the inhomogeneous ODE system 

(6.4) X{t) = A{t)X{t) + F{t) 
with initial condition X{0) = is given by Duhamel's formula 

(6.5) X{t) = [ S{t,t')F{t')dt' 
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Proof of Lemma \6.1\ Let a 



where g = g{t) is given by 
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- a and v = v — v; these perturbative functions satisfy 




5(0) 
^(0) 








/(«) - /(a) 



g 



a — a 

which is (in particular, uniformly bounded). Set 



if a 7^ a 
if a = a 



A{t) 



1 

qg{t) 



a{t) 
v{t) 



in (16.41) . and appeal to Duhamel's formula (16.51) to obtain 



(6.6) 



d{t) 

m 



dt' 



Apply the Gronwall estimate (16.31) with 

w{t)_ 

to conclude that 



S{t + t\t') 



"el(t')" 




'ho 






/2(t'). 
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el(t') 



< V2 



-i/2exp(|g|i/2(t-f)y 
exp(|g|i/2(t-t')) 



max(|g|V2|ei(t')|, 1^2(01) 



Feed this into (16.61) to obtain that on [0,T] 



i{t)\ < V2T 



exp(|g|i/2T) 



1/2 



sup max(|g|^/^|ei(s)|, |e2(s)|) 

0<s<T 



\vit)\<V2Texp{\q\^/^T) sup max{\q\^/^\ei{s)\,\e2{s)\) 

0<s<T 

Taking T < 5|g|~^/^ log(l/|g|), we obtain the claimed bounds. 

7. Main theorem and proof 



□ 



Here we put all the components together and give a stronger version of Theorem [TJ The 
basic procedure is the iteration of Lemmas 15.31 and 16.11 which can roughly be described as 



follows: if the conditions (15.51) hold, and the initial condition satisfies < , say, 

then on the interval ||i/;(t)||j:/i < 2|g|-'^~'^. That means that the evolution of the 

parameters g(t) ^ G is close to the evolution using the effective Hamiltonian, in the way 
that makes Lemma [6.11 applicable . But that gives us a lower bound on tj+i for which (15.51) 
holds on closing the bootstrap loop. 

More precisely, we have 
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Theorem 2. Suppose \q\ <^ 1 and \vq\ ^ 1. Let u solve 

idtu + dlu — q6o{x)u + \u\'^u = 

with initial data Uo{x) satisfying 

\\uo-e'""vi*-ao)\\m<C\q\. 

Then, for times < t < 6{vq + log(l/|g|), the smallness condition in Lemma \3. 1\ 

is met, and thus there are parameters fi, v, 7, a satisfying the symplectic orthogonality 
conditions stated there. Furthermore, we have 

\\u - /ie'"V^r/(/i(x - a))||Hi < c\q\^-'^^ 
Moreover, if a, v, 7 solve the ODE system 

(7.1) t = v, V = ~qd^rj^{a) , ^ = ]^v'^ + ]^- qri^{a) + ]^qad^ri^{a) . 

with initial data (ao,t'o,0), then 

\a - a\ < c\q\^~^\ I7 - 7I + |t; - t;| < c|g|^"^^ |/i - 1| < c|gp"^ 

Proof. The equations (17.11) imply the conservation of energy 

11 11 

-V + -qri (a) = -v^ + -qrj (aq) 

from which we obtain the bound 



(7.2) \v\ < ^Jvl + 2\q\. 

Let 

e{tf= sup \\w{s)\\l, 

0<s<t 



By Proposition 15.21 

(7.3) \a-v\ + \v + ^qd^T]'^{a)\ < Coi\q\\\w\\ + \\w\\hi + ||w||f^i) . 

Let ti with T > ti > be the maximal time for which the assumptions of Lemma 15.31 
(15. 5p hold with i = 0. Then by Proposition 15.41 with n = 1, we have e^(ti) < lo'P"''- The 
estimate (17.31) implies (16. ip in Lemma [6?T] for t G [0, ti], with /(a) = —dxf]{a)/2. By Lemma 
16. II and (17. 2p . we have 



max \v(s)\ < 2\ Vn + 2\q\ . 

0<s<ti \ " 

Reviewing (15.51) . we now see that 

T>ti >C4(^;o' + 2|g|)-^/^ 
where C4 depends only on the implicit absolute constant in (15. 5p . 



30 J. HOLMER AND M. ZWORSKI 

Now let ^2 with T > ^2 > be the maximum time such that (15. 5p holds with i = 1. 
Then by Proposition 15.41 with n = 2, we have 



e\h) < \q 



2-5 



By (17.31) . we have that (16.11) in Lemma [6.11 holds on [0,^2]- By Lemma [6.11 and (17. 2p . we 
have 



max \v(s)\ < 2\/vn + 2\q\ . 

0<s<i2 V 

Reviewing ( 15.51) . we now see that 

1^2 -til > c^{vl + 2\q\)-^/\ 

with the same C4 as in the previous paragraph. 
Continue until the nth step is reached, where 

n = — 1 , 

logc* 

which is the most allowed in Proposition 15. 4[ But now we know that 

T>t^>c6iv', + 2\q\)''/Hogil/\q\), 

and that on [0, 

\a-a\< \q\^~^^\og{l/\q\), \v - v\ < |g| f-^"^ log(l/|g|) 
We also have from Proposition 15.21 
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< 2\\w\\'jji + \q\ IIj/i 



Subtracting the equations for 7 and 7 and using that ||w|| < |gp ^, we obtain 

I7 — 7I < — v"^] + |g||?7^(a) — ?7^(a)| + |g||a — a|?7^(a) + |g| |a| |9^?7(a) — dlr]{a)\ 

< {iq^qf^-'' + \q\\q\'~'' + \q\\q\'-'') log(l/|g|) + Ig^"^^ log^(l/|g|) 
<|gp-2^1og(l/|g|) 

Since we restrict to times t < ^Igl^"*^/^ log(l/|g|), we integrate to obtain I7— 7I < |g|t~^''. □ 



Remark. There remains the case of initial velocities, wq, which are not small. When 
|g| ^0 and fo > is fixed, the dynamics is not interesting and the solution can be 
approximated by the solution with q = 0, that is by the propagating soliton ( 11.101) . The 
proof of that follows from the arguments of [HI §3.1]. 
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